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Закономерности поведения оценок точности 
полиномиального преобразования 
геометрической коррекции изображений 
космической съемки 


На основе анализа основных соотношений полиномиального преобразования изображений космической 
съемки наземных объектов, используемого в процессе их топографической привязки, показана возможность 
оценок ее точности. Предложена методология соответствующих расчетов и моделирования, позволяющая 
как выполнять расчеты в конкретных задачах, так и вскрывать пути снижения ошибок, а также получить 
потенциальные характеристики точности, сняв имеющие место ограничения в распространении достигну- 
тых результатов на любой практически возможный случай. 


Введение 


Тематическое использование изображений космической съемки, как правило, 
имеет целью выяснение положений тех или иных природных или антропогенных 
образований на поверхности Земли. Для каждого из таких применений допустимые 
отклонения соответствующих координат от истинных значений могут существенно 
отличаться. При этом очевидно, что существует целый ряд специфических областей 
использования результатов дистанционного зондирования Земли, где востребованы 
наиболее высокие точности знания топографического положения объекта. К ним, на- 
пример, могут быть отнесены задачи предотвращения или устранения последствий 
чрезвычайных ситуаций, практика так называемого «точного земледелия», а также 
такие наиболее критичные к ошибкам топопривязки проблемы, как сопоставление 
изображений разновременной съемки. 

В целом возникновение ошибок привязки изображений спутниковой съемки может 
быть обусловлено: 

— особенностями применяемой аппаратуры съемки земной поверхности; 

— способом устранения отличий в пропорциях изображения дистанционной съемки 
и соответствующего участка земной поверхности. 

Ограничения здесь связаны с невозможностью безошибочных знаний параметров 
внешнего и внутреннего ориентирования аппаратуры в момент съемки. В то же время 
моделирование полета носителя этой аппаратуры (космического аппарата) и условий 
съемки с некоторым объемом сопутствующих измерений дает возможность получать 
координаты привязки с достаточно низкими ошибками [1]. 

Согласно имеющимся упоминаниям [2-5] «минимальные уровни обработки» кос- 
мических снимков (например, Цопо$ Рго и ОшсКкВпа ОпогесННеа), предлагаемые ком- 
паниями-поставщиками, предполагают достижение среднеквадратической ошибки (СКО) 
их привязки не хуже 5 — 7 м. Однако пользователю придется заплатить за такие снимки 
цену, превышающую стоимость «сырых» изображений в 2 - 3 раза. Если же вести речь 
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о снимках, имеющих СКО привязки до 2 — 4 м (уровни обработки Ргеслуюп, Ргес151оп 
убегео), то соответствующее увеличение достигает уже 6 раз. Причем в ряде случаев для 
обеспечения указанных точностей ортотрансформированного изображения компанией- 
поставщиком от пользователя требуется предварительное предоставление реперных 
(опорных) точек, что наряду с высокой стоимостью еще и существенно удлиняет сроки 
получения материалов съемки [3]. 

Поэтому рассмотренные особенности позволяют считать [3], что пользователи 
могут заказать изображения первого уровня (в 5,5 раза дешевле продукта Ргес1$1оп), 
а дальнейшую обработку произвести самостоятельно. При этом наиболее распростра- 
ненным в тематической обработке изображений космической съемки является при- 
менение полиномиальной геометрической коррекции. Возможности ее достаточно 
широки в выборе параметров осуществления, что в конечном итоге отражается на 
степени повышения интересующей точности плоского геометрического сведения со- 
поставляемых изображений. 

В то же время, приступая к такому виду обработки, исполнитель должен иметь не 
просто представления о характере зависимости ошибок привязки изображений косми- 
ческой съемки от тех или иных факторов, а знания и реализуемые возможности оценок 
и контроля точностей пространственного преобразования, достигаемых в процессе вы- 
полнения этого этапа геометрической коррекции. На создание соответствующих условий 
и ориентирована настоящая статья. 


|1 Аналитические возможности оценки ошибок 
полиномиальной геометрической коррекции 


В основе полиномиальной геометрической коррекции лежит очевидное преобразо- 
вание растровых координат изображения в плоскую прямоугольную систему топографи- 
ческих координат: 


х, =ОГ.У; 
5 сх’т 
ЯУ=У\У у, (1) 
гдеу=| 1х а а р у’ "_ векто построенный 
д и р рр Ар о Ур р Льр Ур: р, тр 
на растровых координатах пикселя изображения, пересчет которых в топографиче- 
ские координаты указанным способом осуществляется (зависит от порядка регрессии); 
х,, у, - растровые координаты пикселя, топографические координаты ко- 
торого х,, у, являются конечной целью преобразования; 
о Е р а хо ^ — 
О = [4 4..2. | М=|® я = статистические (в общем случае) оценки 
коэффициентов полиномиального преобразования. 

Следует сразу обратить внимание, что компоненты векторов У представляют 
собой детерминированные значения. В то же время оценки векторов преобразований 
О, \ ‚ в основе определения которых лежат регрессионные зависимости между раст- 
ровыми и топографическими координатами реперных точек, определенных на местности 
с ошибками, являются случайными величинами. Следствие этого — ошибки оценок 
прямоугольных координат каждого из пикселей изображения х,,у,. 

Отсюда согласно (1) дисперсия возникающих ошибок для координаты х, (для 


координаты У: аналогична) может быть определена так: 


р.) = маг у) 60" -у] = мб" -^0- 0" у} = у" .м(б- 0 у=у".|К|-5, (2) 
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где ЛО — означает разность случайной величины вектора преобразования О и 
его математического ожидания (центрированный вектор значений); 
|К| — ковариационная матрица коэффициентов полиномиального преобразо- 
вания. 
Присутствующая здесь (2) ковариационная матрица полиномиального преобразова- 
ния согласно [6], [7] выражается следующей зависимостью: 


[К[= о» 6. (3) 


Здесь О’рг — среднеквадратический разброс ошибок (СКО) прямоугольных ко- 


ординат привязки реперных точек в топографической системе; 
При этом сама матрица |С| в данной задаче строится на основе использовании 


матрицы предикторных переменных № ь представляемой аргументами полиномиального 


преобразования: 
ТР 
[<= МГУ (4) 
где |\ = У Уз ро У и 
т 
Ур, = | 1 Х ры Уря р и Х ры ь Уря т ИА п . Ур Х ры у У ..:] — компо- 


нента вектора, построенная на полиномиально представленных предикторах растро- 
вых координат 1-й реперной точки (1= 1, т ). 


В итоге, с учетом (3), дисперсия координат любого пикселя преобразуемого 

изображения (2) принимает окончательную запись: 
В, )= оу" |<" у. (5) 

Таким образом, последнее выражение позволяет уже более конкретно судить о 
том, что дисперсии оценок ошибок топопривязки координат геометрически преобразо- 
ванного изображения О(х,), 2(у, ) зависят явно от дисперсии реперных точек (02) и 
неявно 

— от порядка полиномиального преобразования, соответствующего степени ис- 
кажений изображения при съемке (через растровые координаты расставленных 
реперных точек); 

— от числа реперных точек, формирующих матрицу предикторных переменных 
[У] и через нее — матрицу |С 


2 


— от значений растровых координат, определяющих матрицу |С 


‚ а следовательно, 


от расстановки реперных точек. 
Отсюда очевидна возможность получения соответствующих количественных 
данных при любых сочетаниях параметров явных и неявных связей. 


2 Представления о системе исходных данных, 
действиях, моделировании и решаемых задачах 


Создание представлений о величинах ошибок и определяющих их факторах 
обеспечивает построение системы задания данных, взаимосвязанных расчетов и мо- 
делирования, основанное на варьировании условий и параметров реализации выра- 
жения (5). Структурная схема такой системы может быть представлена, как показано 
на рис. 1. 


2 


1 НОЯ Е 
Здесь без ущерба для общности считаем о. =0‘ьр брт: 


«Штучний 1нтелект» 32011 209 


Белозерский Л.А., Мурашко Н.И. 
Задание (выбор) числа 
ренерных точек Иру 


Задание размеров изображения Расстановка реперных Формирование ошибок 
и пространственной точек (определение расстановки реперных 
разрешающей способности растровых координат при точек (датчик случайных 
съемки расстановке) чисел) 


Формирование, расчет и 
обращение матрицы |®] 


Задание вектора растровых Задание СКО топопривязки 
координат У преобразуемого Расчет СКО ошибок 


изображения (выбор из числа й РЕ.) м [26;,} 


возможных) 


реперных точек О’рг 


Графическое оформление 
заданий и расчетов 


Рисунок 1- Структурная схема системы задания, расчетов и моделирования 
в оценке ошибок топопривязки пространственно преобразованного изображения 


Рассмотренная структурная схема определяет полный перечень исходных данных 
в комбинации с взаимосвязанными действиями, обеспечивающими оценку точности 
привязки любой точки реального изображения космической съемки. Получаемые здесь 
ошибки — это результат полиномиальной пространственного преобразования для за- 
данной расстановки реперных (опорных) точек, их числа и точностей топопривязки. 
Причем, располагая точно известными характеристиками перечисленных данных, модель 
оказывается полностью идентифицированной для достижения указанной цели. Соответ- 
ствие натуре при этом обеспечивает, как использование одних и тех же данных не- 
посредственно для геометрической коррекции и ее моделирования, так и применение 
единой для них аналитической основы в преобразованиях координат и в определении 
ошибок. Естественно, что указанное соответствие (в том числе и при статистическом 
задании ошибок) сохраняется в случаях необходимости решения задач не только кон- 
кретного расчета, но и определения путей улучшения результата и поиска закономер- 
ностей при изменении: 

— размеров изображения; 

— пространственной разрешающей способности съемки; 

— растровых координат расстановки реперных точек; 

— точностных характеристик топографического позиционирования реперных точек; 

— случайных ошибок расстановки реперных точек при их интерактивном раз- 
мещении на изображении относительно мест точного положения на местности. 

То есть, изменяя в интересующих пределах перечисленные параметры, можно 
оценивать реакцию на них ошибок привязки геометрически корректируемого изобра- 
жения космической съемки во всем факторном пространстве внешних воздействий, 
определяемых вышеприведенным перечнем. Остается понять, насколько индивидуаль- 
ны получаемые реакции в каждой конкретной комбинации приведенных параметров, 
что и позволяют выполнить имитационное моделирование любой практически име- 
ющей место ситуации согласно приведенной схеме (рис. 1). 
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3 Расчет параметров, его ограничения 
и общие закономерности 


Оценка ошибок геометрической коррекции изображения космической съемки тре- 
бует некоторой начальной определенности параметров. Для конкретизации последующего 
анализа поведения определяемых характеристик используем следующий пример: 

— размеры изображения Мс х М, = 2000 х 2000 пикселей; 


— пространственная разрешающая способность съемки 5 м; 

— расстановка реперных точек — геометрически равномерная; 

— ошибки расстановки отсутствуют; 

— Порядок регрессии топографических координат по предикторам растровых 
координат. 

На рис. 2в условиях указанных ограничений приведено семейство зависимостей 


СКО (с„,=./Р(,) или о, =./(у,) ) топографической привязки северо-западного угла 
геометрически откорректированного изображения х„,У„ для первого и второго порядков 
регрессии. При этом количественно отличающиеся наборы реперных точек (Ир = 5, 10, 
15, 20, 25, 30 _) и среднеквадратические точности позиционирования каждой из них (1,2 и 
5 м) выбраны произвольно". 


=. СКО 
ы =5 ско = р еперной 
е 2 д реперной вес 4 р а 
58 №. точки = 1 м ЗЕз А, 
т 3 оф 
55 о ех.2 
5 52 то = -.-.-ско И Бе И м 
©) © | Нл ГР реперной оз 1 реперие 
бо точки = 2 м о=0 токи 
1234656 т 2 3 5 
_--- ско —--- СКО 
Число реперных  реперной Число реперных реперной 
точек/5 ОИ точек/5 - 1 [м] точки = 5 м 
а) 6) 


Рисунок 2 — СКО топопривязки преобразованного изображения 
для отличающегося числа и точностей определения координат реперных точек: 
а) для линейной регрессии, 6) для квадратичной регрессии 


Приведенные результаты, прежде всего, демонстрируют, что для линейной ре- 
грессии среднеквадратический разброс рассматриваемой ошибки снижается практически 
вдвое в принятом диапазоне изменений числа реперных точек при любой из точностей их 
топографического позиционирования. В то же время для регрессии второго порядка 
такого влияния числа реперных точек на ошибку привязки изображения не наблюдается 
при более высоком ее относительном уровне. 

На первый взгляд, полученные оценки кажутся ограниченными, характерными, 
прежде всего, для приведенного случая в силу единственности размера изображения, 
заданного для моделирования. 


1 у 
Произвольный характер выбора объема первоначальных оценок обусловливался необходимостью наблюдения 
имеющих место закономерностей при всегда имеющей здесь место возможности дополнить состав реперных 
точек и задать интересующие дополнительно ошибки их топопривязки на местности . 
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Однако вопреки этому, повторяя расчеты, но при увеличении изображения (с со- 
ответствующим изменением координат реперных точек), например, в 10 раз (то есть до 
20000 х 20000 пикселей), достаточно просто обнаружить, что точности привязки 
остаются теми же, что и приведены на рис. 2. Причем к такому же результату приходим 
и с уменьшением изображения. Одновременно с равным успехом к нему же ведут отли- 
чающиеся изменения размеров его по каждой из координат. 

Объяснение получаемого результата следует из пропорциональности определи- 
теля (2) матрицы |<] (4) и ее алгебраического дополнения 4, размеру изображения. 


В результате компонента обращенной матрицы |С| не изменяет своего значения: 


си =А,/Р=сопзй. 

При этом вектор растровых координат точки привязки изображения (северо-за- 
падный угол изображения) выглядит для линейной и квадратичной регрессии доста- 
точно просто: 


УЕ О 0] и У=р ООО (]. 


Соответственно среднеквадратическая ошибка топопривязки согласно (5) ока- 
жется одинаковой для любого увеличения изображения И равной: 


Ст, — Орг уу" СГ и [6 ° 


Наряду с установленной особенностью преобразования координат легко убедить- 
ся, что увеличение или уменьшение пространственной разрешающей способности съемки 
при постоянном размере изображения не оказывает никакого влияния на относительное 
расположение расставленных реперных точек, а значит, и на проведенные оценки. Эти 
изменения в анализируемом преобразовании оказываются равносильны изменениям раз- 
меров изображения. 

В результате можно считать, что первое и второе из введенных начальных огра- 
ничений моделирования, связанные с выбором геометрических параметров изображения 
и разрешающей способности, могут быть сняты, а полученные зависимости точностей 
привязки пространственно преобразуемого изображения (рис. 2) следует относить к не- 
изменным характеристикам, потенциально достигаемым (оценка снизу) для любого по 
размерам изображения, для любой разрешающей способности. Пропорционально увели- 
чивая или уменьшая изображение (до реально полученных размеров в космической 


1 
съемке), в любом практически имеющем место случае заданных величин Орг При со0- 


ответствующем изменении координат идеально расставленных реперных точек приходим 
к тем же зависимостям потенциальной точности привязки. 

Несмотря на такие преимущества, рассмотренные характеристики все-таки ка- 
саются оценки только нижней границы ошибок. Само появление приведенного выше 
определения «потенциальная точность» обусловлено здесь тем, что «чистая» равно- 
мерность расстановки реперных точек может быть обеспечена только при моделиро- 
вании. Практическое достижение такого же качества подготовки геометрической кор- 
рекции реальных изображений космической съемки оказывается невозможным. Здесь 
при стремлении к желаемой равномерности начинает играть роль сложность выбора 
достаточного числа характерных мест расстановки реперных точек, не говоря уже о 
невозможности топопривязки их в каждом из заранее выбранных мест. 


`В силу незначительности принятого здесь для моделирования диапазона изменений СКО топопривязки 
реперных точек возможно также ошибочное отнесение результатов, полученных на рис. 2, к трудностям рас- 
пространения их на любой практически имеющий место случай противоречит выражению (5). Пропорционально 
увеличивая или уменьшая приведенные значений с’ = /р,) в число раз, в которое отличается реально 


интересующее значение Орг от моделируемых, легко удовлетворить любые потребности. 
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Являясь ориентиром при организации обработки реальных изображений косми- 
ческой съемки и пределом в их привязке, потенциальная точность не исключает оценку 
реально достигаемых значений ошибок. Для практически имеющих место расстановок 
реперных точек ничто не препятствует их определению на основе использования той 
же, уже рассмотренной методологии. 

Таким образом, сняв предшествующие ограничения, пока остается необходимым 
считаться с очередным из перечисленных в начале раздела ограничений: моделирование 
выполнялось в предположении, что ошибки расстановки реперных точек отсутствуют. 

Насколько серьезно такое допущение, позволяет установить теперь уже статисти- 
ческое моделирование, предусмотренное согласно структурной схеме рис. 1. Для этого 
достаточно задействовать датчик случайной расстановки реперных точек оператором. 
Причем в качестве закона распределения соответствующих ошибок оператора можно 
принять равномерное распределение в диапазоне его способностей. Выбирая для исследо- 
ваний два крайне возможных диапазона этих ошибок (+100 пикселей — недостаточно 
внимательный оператор и +3 пикселя — аккуратный оператор), удается проследить ин- 
тересующие закономерности. 

Согласно изложенным условиям на рис. 3 приведена статистически получен- 
ная выборка среднеквадратических ошибок (СКО) привязки изображения космиче- 
ской съемки при расстановке 10 реперных точек с точностью знания их положения на 


местности Орг =1м. 
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Рисунок 3 — Реализации СКО привязки изображения при растровой 
расстановке 10 реперных точек с ошибками: а) +100 пикселей; 6) +3 пикселя) 


Качественно оценивая один и другой результат, можно отметить, что относительно 
ошибки привязки изображения 0,667, достигаемой при точной расстановке 10 реперных 
точек с о», =1м (рис. 2а), разброс, обусловливаемый неточностями операторов в их пра- 


ктической работе по ручной расстановке, незначителен. 

В табл. 1 для тех же исходных данных по разрешающей способности и точности 
топографического позиционирования реперных точек сведены статистические оценки 
привязки геометрически преобразованного изображения в диапазонах ошибок расстанов- 
ки, более приемлемых для практики. 

Как можно увидеть по приведенным данным (рис. 2, табл. 1), неточность растро- 
вой расстановки заданного числа реперных точек с ростом ошибок оператора растет не 
существенно. Отклонения от математического ожидания при возрастании диапазона 
ошибок расстановки от +3 до +30 пикселей составляют 0,2% — 2%. 
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Таблица 1 
Диапазон Интервал 
ошибок Математическое Среднеквадратический практически 
расстановки ожидание СКО разброс СКО привязки возможных 
реперных привязки [м] значений 
точек [м] точности 
[пикселей] привязки 
(по выборке) 
[м] 
—30 - +30 0,667 0,0044 0,656 = 0,675 
—20 = +20 0,667 0,0028 0,661 0,672 
—10 = +10 0,667 0,0013 0,663 -0,668 
—3++3 0,667 0,0004 0,665 = 0,667 
Здесь же следует обратить внимание на то, что, как показывает моделирование, 
изменения разрешающей способности для того же диапазона операторских ошибок 


расстановки реперных точек не сказывается на ошибках топопривязки изображения 
космической съемки. 

Таким образом, и начальное ограничение в рассмотренном моделировании, со- 
стоящее в отказе от учета ошибок расстановки реперных точек, оказывается допустимым, 
если не предполагать их необоснованно большими. Например, имея дело со случаем, 
когда ширина проезжей части пересекающихся дорог составляет 8 м, при желании уста- 
новить реперную точку в центре их перекрестка, хорошо отображаемого космическим 
снимком с пространственной разрешающей способностью 1 м, ошибка установки практи- 
чески никогда не превысит 2 — 3 пикселя. Но даже необоснованно выставляя эту точку 
на краю одной из дорог, хотя и в районе пересечения вместо естественного стремления к 
центру, придется иметь дело с максимально возможной ошибкой всего в 4 пикселя. При 
равной тщательности расстановки других реперных точек это не приведет к ухудшению 
ошибки привязки геометрически скорректированного изображения (рис. 2) более, чем 
на 0,2%. 

Аналогично обстоит дело и при ухудшении разрешающей способности. Если в 
отличие от предшествующего случая интересующее пересечение дорог выглядит всего 
лишь как пересечение линий, то в качестве места установки реперной точки не может 
быть принят другой пиксель, кроме того, в пределы которого это пересечение попадает. 


4 Точность геометрической коррекции и степень 
равномерности расстановки реперных точек 


Возвращаясь к вопросу равномерности расстановки реперных точек, придется под- 
черкнуть, что в качестве действий по преодолению ограниченности моделирования в 
соответствующем случае предлагалось выполнение отдельного расчета согласно рас- 
смотренной модели для реально имеющих место условий геометрической коррекции. 
При этом, как следует из рис. 2, точность топографической привязки изображений может 
быть повышена увеличением числа реперных точек и снижением ошибок определения 
их положений на местности. Такой путь, как можно было понять, часто более сложен, 
чем стремление к равномерной расстановке. 

Здесь приходится обратить внимание, что рассмотренная методология подготовки 
данных, расчетов и моделирования наряду с перечисленными способами приближения 
к желаемой точности геометрической коррекции открывает одновременно и возмож- 
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ности для изменения расстановки реперных точек. К открытию ее приводят попытки 
варьирования расположением реперных точек в пределах изображения, а соответст- 
венно на отображаемой части земной поверхности, с оценками точности линейной 
регрессии. 

В процессе исследования различных вариантов моделирования может оказаться, 
что и реально неравномерное расположение реперных точек в некоторых случаях обеспе- 
чивает ту же потенциальную точность привязки северо-западного угла изображения. 
В частности этому соответствует сосредоточение реперных точек в непосредственной 
близости от указанного угла, в то время как в остальной части изображения они вообще 
не расставлялись. То есть графики зависимостей о, (срт,п»ь) В таком случае остаются 


теми же, что и приведены на рис. 2. 

Для расширения представлений о таких ситуациях следует не ограничивать инте- 
ресы оценки точности топографической привязки единственной рассмотренной точкой 
изображения (Х„,У„), играющей важную роль в его применении, а распространить 


расчет на все пространство изображения. В результате можно обнаружить закономерные 
изменения точностей привязки изображения космической съемки по всей его площади. 
Так, на рис. 4 приведены ошибки топографической привязки сетки пикселей изображе- 
ния, полученные моделированием линейной регрессии для 10 равномерно расставленных 
реперных точек и орг =1м. 


Получаемая поверхность ошибок в физической аналогии представляет собой 
«платок, подвешенный закреплением его на равных высотах за четыре угла». Причем 
вертикально вниз от каждой точки закрепления располагаются углы рассматриваемого 
изображения, занимающего горизонтальную часть рисунка. 

Здесь можно увидеть, что максимальная точность привязки, обеспечиваемая гео- 
метрической коррекцией, достигается в точке максимального провисания «платка», сов- 
падающей с центром изображения (0,33 м), в то время как значение 0,67 м, приведенное 
на рис. 2 для одного указанного выше угла, характерно и для всех его угловых положений 
(«высоты закрепления»). 

Причем, увеличивая число реперных точек, принимающих участие в такой оцен- 
ке, придется увидеть, что представленная на рис. 4 поверхность перемещается вниз по 
оси 7, (уменьшаются высоты точек закрепления «платка» и точки его максимального 
провисания, соответствующие ошибкам топопривязки). 

Тем самым отражается повышение точностей для каждого из угловых положений 
изображения и, в частности, центра его. 
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Рисунок 4 — Среднеквадратические ошибки топопривязки сетки пикселей 
изображения при равномерной расстановке 10 реперных точек (©; = 1м) 
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Закономерность такого поведения оценок на сетке пикселей легко проверяется 
теоретически, если попытаться определить экстремум среднеквадратической ошибки 
топографических координат как функции растровых координат, представляемой вы- 
ражениями (2) или (4), в плоскости изображения: 


= 2550 
ОХр р 


Тогда с учетом особенностей обращенных матриц |С| для первого порядка ре- 


грессии в соответствии с приведенными экспериментальными данными получим, что 
экстремум соответствует минимуму СКО и достигается в точке с растровыми ко- 
ординатами, определяемыми через коэффициенты корреляции, как составляющие ко- 
вариационной матрицы |К| (2): 
шт ы В 12. шт . Е К 
х = агв т Б(х,) Е У, = агвтит Б(%,) = 
р 22 ыы 33 
При конкретных значениях коэффициентов корреляции линейной регрессии, 
получаемых в процессе моделирования, и равномерной расстановке реперных точек 
используемого изображения 10000 х 10000 м получим: 


в =-5000 м; у = 5000 м. 


Теперь можно воспользоваться моделированием для объяснения упомянутой выше 
ситуации сохранения потенциальной точности с Хоз Оуюо ПРИ группировании реперных 


точек в пределах некоторой части изображения, примыкающей к северо-западному углу 
его. В частности, для этого при моделировании достаточно уменьшить, например, в 2 раза 
значения заданных координат реперных точек (равномерного размещения) при сохране- 
нии размеров моделируемого изображения. 

В результате такого действия область расстановки реперных точек окажется за- 
нимающей левую верхнюю часть пространства изображения, составляя четверть его, 
в пределах которой те же 10 реперных точек размещаются по-прежнему равномерно. 
Соответственно этому на рис. 5 приведены ошибки топографической привязки сетки 
пикселей всего изображения. 


СКО топопривязки 
пикселей 
изображения [м] 


Координаты по оси 
У изображения 
Ру=у*2000 [м] 


Координаты по оси Х 
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Рисунок 5 — Среднеквадратические ошибки топопривязки сетки пикселей 
при сосредоточенной равномерной расстановке 10 реперных точек 


в левой верхней четверти изображения (©: =1м) 
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Приведенные результаты в отличие от предшествующего случая позволяют обнару- 
жить, что «платок» оказывается подвешенным за углы на отличающихся высотах. То есть 
— среднеквадратическая ошибка топопривязки северо-западного угла изображения 
(левый передний угол поверхности на рис. 5) сохранила свое значение (с„,= 0,67 м), 


соответствующее равномерной расстановке реперных точек в пределах всего изобра- 
жения (рис. 4); 

— среднеквадратические ошибки других угловых положений изображения су- 
щественно возросли и отличаются между собой; 

— минимальная ошибка топопривязки сохранила свое значение (0,33 м); 

— местоположению минимальной ошибки соответствуют растровые координаты 


изображения х"" = 2500 м; ум" = 2500 м. 


Более тщательное изучение рассматриваемого результата дает возможность увидеть: 

— минимальная ошибка топопривязки (0,33 м) соответствует не центру всего изо- 
бражения, а центру области сосредоточения равномерно распределенных реперных точек 
(5000 х 5000 м); 

— на углах полученной здесь области сосредоточения реперных точек ошибки 
топографической привязки пикселей равны между собой и их СКО составляет 0,67 м, 
как это было для всего изображения при равномерной расстановке реперных точек в 
его пределах. 

На основе полученных результатов придется вспомнить, что полиномиальная ап- 
проксимация зависимости между растровыми и топографическими координатами для 
регрессии первого порядка строится на данных области расстановки реперных точек. 
Полиномиальное определение координат каждой точки изображения по восстановленной 
указанным образом регрессии представляет собой интерполяцию. При ее выполнении 
достигается минимум ошибки в центре этой области и равные значения в равноудален- 
ных и симметричных относительно центра точках. Чем более удалены такие точки, тем 
больше СКО их привязки. 

Центр, определяемый по минимуму ошибки интерполяции, всегда соответствует 
центру области равномерной расстановки. 

В случае, если область равномерной расстановки реперных точек по занимаемой 
площади меньше площади изображения, то для определения ошибок привязки точек 
изображения, лежащих вне указанной области, полиномиальная интерполяция пере- 
ходит в экстраполяцию с естественным ростом ошибок. 

Таким образом, с какой бы расстановкой реперных точек ни пришлось иметь дело, 
определение соответствующего центра (по минимуму ошибки привязки сетки пикселей 
изображения) и его смещения относительно истинного центра интересующего изобра- 
жения для регрессии первого порядка позволяет: 

— установить степень равномерности имеющей место расстановки в пределах 
всего изображения; 

— дополнить имеющийся состав реперных точек для достижения большей равно- 
мерности; 

— сопоставляя ошибки привязки углов изображения с потенциально возможными, 
определить степень удаления области расстановки от соответствующих углов изобра- 
жения. 


Заключение 


Предпринятый в статье поиск закономерностей поведения ошибок полино- 
миальной коррекции изображения космической съемки основывался на расчетах и 
моделировании таких оценок во всем факторном пространстве изменений соответ- 
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ствующих внешних воздействий. Предложенная методика оценки ошибок полиномиаль- 
ного преобразования координат открыла возможности углубленного анализа причин 
их возникновения, наглядно и количественно демонстрируя их связи с такими факто- 
рами, как 

— число реперных точек, участвующих в восстановлении регрессионной зави- 
симости между растровой и топографической системами координат; 

— точности привязки координат реперных точек в топографической системе; 

— порядок регрессии связи топографических координат точек изображения с их 
растровыми координатами; 

— степень равномерности распределения реперных точек на изображении, соответ- 
ствующем отображаемому участку земной поверхности. 

В процессе исследований, результаты которых приводятся в статье, для каждого 
порядка регрессии установлена независимость потенциальной точности топографи- 
ческой привязки координат изображений космической съемки от 

— размеров изображения; 

— разрешающей способности космической съемки; 

— ошибок растровой расстановки реперных точек на изображении космической 
съемки. 

На основе исследования результатов моделирования предлагается методология, 
позволяющая пользователю в процессе тематической обработки данных космической 
съемки: 

— располагать заранее детерминированными графическими или табличными дан- 
ными для оценки возможного нижнего предела ошибок (потенциальная точность) при- 
вязки полиномиально преобразуемого изображения применительно к его параметрам; 

— получать оценки ошибок топопривязки геометрически трансформированных 
изображений космической съемки для реального числа, заданной расстановки репер- 
ных точек и знаний ошибок их положений на местности; 

— определять требуемые ошибки позиционирования на местности и количество 
используемых реперных точек при необходимости повышения точности привязки 
изображений космической съемки; 

— оценивать равномерность расстановки реперных точек и при необходимости 
повышения точности привязки изображения корректировать ее на основе оценки 
ошибок топографической привязки сетки пикселей изображения в целом. 
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Л.А. Ылозерський, М1. Мурашко 

Закономфнеть поведнки ощнок точност! полном1ального перетворення геометрично! корекцй 
зображень косм!чно! зйомки 

На основ! аналзу основних сшведношень полномпального перетворення зображень косм!чно! зйомки 
наземних об’екив, використовуваних в процес! 1х топографлчно! прив’язки, показана можлив1сть ощнок й 
точност!. Запропонована методоломя вдповдних розрахунюв 1 моделювання, шо дозволяе як виконувати 
розрахунки в конкретних завданнях, так 1 розкривати шляхи зниження помилок, а також отримувати потен- 
щин! характеристики точность, знявши обмеження, що мають мсце, в поширенн! досягнутих результатв на 
будь-який практично можливий випадок. 


Г.А. Веозег К), М. МигазйКо 

ВесшагИу о! Вевауюг 01 Ассигасу ЕзИтайоп аё Роупопи а] Тгап$огтайоп о{ Ппабе Зрасе Зигуеу 
Сеосгарше Согг@айоп 

Оп Фе Баз15 оРапа]уз1$ оЁ Фе Базс ргорогйоп оЁ ро]упопта] гапзЮптайоп оЁ птазез оЁ зрасе зигуеу оЁ Фе 1ап4 
обес, \ЫсЬ 15 изе4 фамие Фе ргосез$ оЁ Фе ‘юроэтарыс теЕтгепсе оР фезе обес, роззОПИу оЁ ассигасу 
езйтланопз 1$ зпо\т. ТЬе теводоюзу оЁ зресла| сасшайоп$ ап то4дейпе Юг ргасйса| сасшайоп, Юг еггог 
тедасноп ап4 а15о Юг сете роепва| свагасет$Нс$ оЁ ассигасу 15 оНегеа. Тиз Фе гезёлсНоп ш ‹ззетттаноп оЁ 
{Фе асшеуеа гези $ ш апу роз е зИлайоп 15 гепоуе4. 
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